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ABSTRACT 


Composite  materials  are  seeing  increased  use  in  structural  applications  because  of 
their  various  benefits.  When  composite  structures  are  employed  in  a  water  environment, 
their  dynamic  responses  are  greatly  affected  by  the  fluid  m  edium.  Water  density  is 
comparable  to  many  composite  materials  and  the  effects  of  fluid-structure  interaction  on 
dynamic  behaviors  of  composite  structures  are  significant.  The  effects  of  fluid- structure 
interaction  include  changes  of  frequency,  magnitude,  en  ergy  dissipation,  etc.,  of 
structural  characteristics.  Hence,  it  is  critical  to  understand  the  fluid-structure  interaction 
of  composite  structures  subjected  to  dynamic  loading  in  water  environments.  This  work 
focuses  on  finding  param  eters  affecting  the  transient  dynam  ic  responses  of  com  posite 
structures.  Coupled  fluid-st  ructure  interaction  analyses  of  com  posite  plates  are 
conducted  numerically,  using  finite  element  models,  including  various  parametric  studies. 
The  results  are  compared  to  those  of  dry  structures  to  identify  the  role  of  each  parameter. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Composites  are  seeing  increas  ed  use  in  m  aritime,  aerospace  and  a  utomotive 
structures  used  in  both  civil  and  m  ilitary  applications.  Early  uses  of  c  omposites  were 
limited  to  secondary  structures;  however,  as  knowledge  and  understanding  of  mechanical 
characteristics  of  composites  has  grown,  more  primary  load-bearing  structures  have  been 
fabricated.  In  recent  years,  la  rge  composite  structures  have  been  incorporated  into  naval 
vessels  to  increase  operational  perf  ormance  wh  ile  lowering  ownership  costs  [1].  For 
example,  carbon-fiber  com  posite  material  pr  ovides  high  strength  and  stiffness  while 
maintaining  low  weight,  which  in  tur  n  translates  to  increas  ed  fuel  econo  my  or  increased 
payload.  A  further  advantage  of  com  posites  over  m  etals  is  lower  m  aintenance  and 
resistance  to  corrosion,  making  composites  very  desirable  for  maritime  applications.  The 
use  of  com  posites  in  engineering  components  has  initiated  num  erous  studies  to  analyze 
structural  components  fabricated  from  various  com  posites  rather  than  traditional  m  etals. 
While  com  posites  provide  advantages  over  m  etals,  they  als  o  com  e  with  com  plex  and 
challenging  engineering  problems  for  analysts  and  designers  [2].  Because  the  stru  ctural 
behavior  is  im  pacted  by  Flui  d  S  tructure  In  teraction  (FSI),  th  is  work  focuses  on  the 
implications  of  utiliz  ing  com  posite  structu  res  in  m  aritime  applications  below  the 
waterline. 

B.  LITERATURE  SURVEY 

It  is  critical  to  assess  the  structu  ral  be  havior  of  com  posite  structures  used  in 
marine  applications  beneath  the  wa  terline  wher  e  FSI  p  lays  an  im  portant  ro  le  on  the 
dynamic  response  and  failure  of  the  submerged  composite  structure.  Because  composite 
structures  are  m  uch  lighter  than  m  etallic  structures,  the  effect  of  FSI  is  m  uch  greater. 
Many  polymer  composite  materials  are  on  ly  a  few  times  heavier  than  water;  the  refore, 
the  added  mass  effect  of  the  fluid  becomes  critical. 

Numerous  studies  hav  e  exam  ined  the  e  ffect  of  FSI  f  or  m  etallic  struc  tures, 
especially  f  or  underwater  explosive  loading  [3]— [10].  Som  e  works  are  experim  ental 
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studies,  while  others  are  num  erical  work.  A  few  studies  exam  ined  FSI  for  com  posite 
structures  subjected  to  underwater  explosion  [1 1]— [17]. 


C.  OBJECTIVES 

This  work  investigates  the  effects  of  the  surrounding  fluid  on  dynam  ic  responses 
of  composite  structures  subjected  to  a  mechanical  loading  via  applied  concentrated  force, 
uniform  pressure  and  i  mpact.  The  research  exam  ines  several  param  eters  affecting 
transient  dynam  ic  responses  of  subm  erged  composite  s  tructures  to  iden  tify  m  ajor 
controlling  param  eters  of  FSI.  Thi  s  resear  eh  focuses  on  com  putational  m  odeling  of 
coupled  fluid-structure  interaction  analys  es  of  com  posite  structures — specifically 
plates — under  water  for  various  param  etric  studies.  Results  are  norm  alized  to  those  of 
completely  dry  structures  to  illustrate  the  role  of  each  parameter  on  FSI. 
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II.  COMPUTATIONAL  MODEL 


A.  MATERI  AL  SPECIFICATIONS 

The  composite  material  used  in  this  study  is  an  e-glass  woven  fabric  with  a  plain 
weave  fiber  architecture  and  vinyl-ester  resin.  The  com  posite  has  elastic  m  odulus  17 
GPa,  Poisson’s  ratio  0.3,  and  density  of  2020  kg/m3.  To  make  a  fair  comparison  between 
dry  and  wet  structures,  any  potential  change  of  composite  material  properties  associated 
with  moisture  absorption  from  water  is  not  considered.  The  steel  used  for  im  pact  study 
has  elastic  modulus  200  GPa,  Poi  sson’s  ratio  0.3,  and  density  8000  kg/m  .  For  a  dry 
structure,  i.e.,  in  a  ir,  there  is  no  spe  cific  modeling  of  the  air  m  edium.  For  m  odels  that 
examine  FSI,  the  water  is  modeled  with  a  density  of  1000  kg/m  ,  and  bulk  modulus  of  2.2 
GPa,  while  water  viscosity  is  neglected. 

B.  FINITE  ELEMENT  MODEL  DEVELOPEMENT 

As  an  initial  step  in  studying  Fluid  Structure  Interaction  effects,  only  linear  elastic 
behavior  is  considered  in  th  is  study.  Solid  m  aterials  are  modeled  using  the  Lagrangian- 
based  finite  element  method,  while  fluid  is  solved  using  the  Eulerian-based  finite  element 
method  [18].  The  com  posite  plate  used  in  the  study  is  thin  (0.002  m  thickness),  having 
an  aspect  ratio  of  at  least  150  (length  to  thickness),  and  nece  ssitates  modeling  through 
shell  elements. 

Due  to  the  thin  com  posite  plate  re  quiring  to  b  e  m  odeled  with  shell  e  lements, 
coupling  between  the  fluid/com  posite  interfaces  presented  a  challenge,  as  th  e  interface 
between  them  needs  to  be  uniquely  defined  by  a  volum  e  or  solid  elem  ents.  To  have  a 
uniquely-defined  volume,  a  stiffened  composite  plate  is  used  to  create  a  unique  volume  of 
composite.  The  stiffened  com  posite  shell  structure  is  composed  of  top  and  bottom  skin 
plates,  coupled  through  vertical  stiffeners.  Each  skin  plate  is  0.3  m  x  0.3  m  and  0.002  m 
thick  and  is  m  odeled  with  e-glass  com  posite.  The  stiffeners  are  m  odeled  of  the  sam  e 
composite  material  with  the  same  thickness  of  the  skin  plates.  Their  sizes  are  0.3  m  long 
and  0.01  m  tall,  and  spaced  every  0.05  m  apart.  The  spacing  between  nodes  of  composite 
model  was  1  cm  ,  such  that  the  0.3  m  skin  pi  ates  have  a  30  by  30  m  esh.  The  stiffened 
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composite  plate  is  depicted  in  Figure  1,  with  the  lines  denoting  the  locations  of  stiffeners. 
With  this  stiffened  composite  plate,  FSI  can  be  investigated  by  comparing  three  different 
cases:  1)  completely  dry,  2)  two-sides  wet,  and  3)  one-side  wet. 


Figure  1.  Stiffened  Composite  Plate  Structure 

Various  parametric  studies,  including  boundary  conditions  and  loading  types,  are 
examined  to  investigate  FSI  effects.  The  edges  of  the  stiffened  composite  plate  use  either 
a  clamped  or  sim  ply-supported  boundary  condi  tion.  The  plate  is  s  ubjected  to  constant 
applied  force  at  the  cen  ter  of  the  top  skin  plate,  equivalent  pressure  loading  over  the 
surface  of  the  top  plate,  or  impact  loading  at  center  of  the  top  plate  from  a  steel  projectile 
at  various  initial  velocities. 

The  Finite  Elem  ent  Models  (FEM  )  were  constructed  in  PATRAN  and  solved 
numerically  using  DYTRAN.  The  com  putations  were  run  using  a  FIPC  cluster  system  . 
The  computational  time  required  to  perform  0.05  second  transient  solutions  varied  from 
approximately  5  minutes  for  the  dry  structure  to  as  much  as  40  hours  for  the  one-side  wet 
structural  model.  The  dry  case  structural  model  has  2,220  elem  ents  and  the  wet  m  odels 
have  up  to  30,000  elem  ents.  The  geometry  used  to  define  the  composite  material  uses  a 
Lagrangian-based  quadrilateral  shape  for  defining  the  shell  elements.  The  geometry  used 
to  define  the  Eulerian-based  fluid  is  composed  of  hexagonal  solid  elements. 
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1.  Dry  Structure 

The  reference  case  throughout  the  study  is  a  completely  dry  structure  using  only 
the  composite  plate  constructed  as  describe  d  previously.  No  sp  ecific  m  odeling  of  air 
surrounding,  and  within  ,  the  void  spaces  of  th  e  stiffened  com  posite  plate  structure  is 
accounted  for,  due  to  its  negligible  effects.  The  dry  structure  dynam  ic  response  is  used 
for  normalization  with  other  cases  to  show  the  effects  of  FSI. 

2.  Two-Sides  Wet  Structure 

A  two-sides  wet  structure  is  used  to  exam  ine  the  influence  of  fluid  (water)  on  the 
response.  It  is  m  odeled  with  the  stiffened  composite  plat  e  embedded  within  a  cube  of 
water.  The  surrounding  fluid  dom  ain  is  much  greater  than  the  composite  plate  structure 
with  a  two  to  one  ratio  of  largest  dimension.  Addition  ally,  the  non-reflectiv  e  boundary 
condition  is  applied  to  the  outside  fluid  boundary.  Although  there  may  be  some  reflected 
waves  from  the  non-reflective  boun  dary  due  to  im  perfect  boundary  condition,  th  e  tim  e 
period  of  interest  for  structural  response  is  too  short  to  include  the  effe  cts  of  reflected 
waves. 


3.  One-Side  Wet  Structure 

A  one-side  wet  structure  is  used  to  simulate  a  condition  in  which  fluid  is  on  one 
side  of  the  plate  while  air  is  on  the  other,  such  as  would  be  encountered  in  construction  of 
a  ship  hull  with  com  posite  p  lates.  To  create  an  air  sp  ace  on  one  side  of  the  stiffened 
plate,  five  additional  rigid  composite  sides  are  added  below  the  stiffened  plate.  The  sides 
are  rigid  sh  ells  composed  of  the  sam  e  composite  m  aterial  and  form  the  volum  e  to  be 
coupled  with  the  surrounding  fluid.  The  one-sid  e  wet  structure  is  depicted  in  Figure  2. 
The  air  volume  between  the  bottom  of  the  stiffened  plate  and  the  bottom  of  the  rigid  box 
is  0.01m  in  height. 


5 


(a)  Box  made  of  a  stiffened  composite  plate  and  five  rigid  sides 


(b)  Composite  box  inside  water 

Figure  2.  Stiffened  Composite  Plate  Supported  by  Rigid  Box  used  in  One-Side  Wet 

Case 
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III.  PARAMETRIC  STUDIES  USING  COMPUTATIONAL  MODEL 


A.  TYPE  OF  BOUNDARY  CONDITION 

Two  different  boundary  conditions  are  applie  d  to  the  stiffened  com  posite  plate, 
clamped  or  sim  ply  supported.  In  reality  it  is  diffi  cult  at  best  to  achieve  a  perfectly 
clamped  boundary  cond  ition,  and  th  us  actual  boundaries  are  a  m  ixture  of  clam  ped  and 
simply  supported.  To  bound  the  dynam  ic  response  of  composite  plate,  both  boundaries 
are  app  lied  individu  ally  to  determ  ine  a  ny  difference  between  FSI  effects.  Any 
experimental  work  done  in  conjunction  with  this  study  will  have  im  perfectly  clamped 
boundaries,  and  thus  the  behavior  will  be  a  mixture  of  both  boundary  co  nditions.  These 
numerical  models  can  be  used  to  understand  the  differences. 

B.  APPLIED  LOADING  TYPE 

The  basis  for  this  study  uses  an  applied  concentrated  force  of  1000N  at  the  center 
of  the  top  skin  plate  to  observe  the  dynam  ic  response  and  detenn  ine  the  FSI. 

Additionally,  an  equivalent  pressure  to  the  concentrated  force  is  also  examined  to  reveal 
any  differences  in  response  from  loading  methods.  Finally,  im  pulse  type  loads  are 
imparted  to  the  com  posite  plate  using  steel  projectiles.  The  steel  projectiles  are  0.3  m 
long,  and  have  either  a  circular  or  square  impact  face  with  area  of  1.6129e'  m  (0.25  in  ). 
The  steel  projectiles  start  2  mm  above  the  top  skin  plate,  and  are  given  an  initial  velocity 
of  1  m/s,  5  m/s  or  10  m/s. 

C.  PLATE  SIZE 

The  basic  stiffened  composite  plate  used  in  this  num  erical  study  consists  of  a 
0.3  m  by  0.3  m  skin  plate.  A  larger  0.5  m  by  0.5  m  skin  plate  model  is  also  examined,  so 
the  differences  in  FSI  can  be  examined  from  increased  spacing  between  supports. 

D.  PLATE  SHAPE 

The  basis  for  this  study  is  the  standard  0.3  m  by  0.3  m  square  stiffened  plate.  To 

examine  the  impact  of  plate  shape  on  the  dynamic  response  and  FSI,  an  equivalent  area 

rectangular  shaped  plate  is  also  modeled  with  dimensions  of  0.2  m  by  0.45  m. 
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E.  COMPOSITE  MATERIAL  PROPERTIES 

Parametric  studies  are  conducted  using  the  basic  0.3  m  by  0.3  m  stiffened  plate  to 
examine  the  effect  of  composite  m  aterial  properties  on  FSI  and  dynam  ic  response.  T  he 
composite  material  is  modeled  with  a  nominal  density  of  2020  kg/m  and  nominal  elastic 
modulus  of  1.7e10  Pa.  Two  different  densities,  approximately  a  50%  reduction  and  100% 
increase  from  the  nom  inal,  are  used  to  inve  stigate  the  change  in  resp  onse;  specifically, 
the  composite  densities  are  1020  kg/m  and  4020  kg/m  .  Two  different  elastic  m  oduli, 
approximately  a  50%  reduction  and  50%  in  crease  from  the  nom  inal,  are  used  to 
investigate  the  change  in  response;  specifi  cally,  the  composite  elastic  m oduli  are  0.7e  10 
Pa  and  2.7e10  Pa. 
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IV.  NUMERICAL  STUDY  RESULTS  AND  DISCUSSION 


A.  METHODOLOGY 

The  basic  methodology  used  to  determine  the  difference  in  dynamic  behavior  of  a 
stiffened  composite  pla  te  is  to  nonnaliz  e  the  tw  o-sides  and  one-side  we  t  cases  with  the 
completely  dry  case.  I  n  this  m  anner,  the  dry  case  is  the  base  re  sponse  and  tends  to 
differentiate  the  particular  changes  due  to  th  e  FSI.  In  general,  th  e  base  case  us  ed  for 
normalization  is  the  completely  dry  plate  with  composite  properties  of  2020  kg/  m3  for 
density,  1.7e  10  Pa  for  elastic  m  odulus  with  clampe  d  edges,  and  a  1000  N  concentrated 
force  applied  at  the  center  of  the  top  skin  p  late.  W  hen  strains  are  exam  ined  the 
nonnalization  is  accomplished  with  respect  to  the  nonnal  x-axis  strain. 

This  method  of  norm  alization  shows  the  tr  ansient  variation  of  various  response 
variables;  such  as  displacement  of  the  central  node  of  top  skin  plate,  strain  energy  and/or 
kinetic  energy  of  the  stiffened  composite  plat  e,  and  stress  or  stra  in  at  one  of  three 
locations  on  the  bottom  skin  plate.  The  numerical  solutions  from  DYTRAN  using  shell 
elements  only  perm  it  stress  to  be  detenn  ined.  Strains  are  calculated  us  ing  the  s  tandard 
stress/strain  transf  onnation  equatio  ns.  In  the  com  putational  m  odel,  stress  (and  hence 
strain  through  transformation  equations)  is  calculated  at  th  e  element  in  center  of  plate  of 
one  quadrant  (this  location  is  termed  ‘center’),  at  an  element  half  way  between  the  center 
and  edge  of  one  quadrant  (this  location  is  te  rmed  ‘side’),  and  at  an  elem  ent  half-way 
between  the  center  and  the  corner  along  a  di  agonal  of  one  quadran  t  (this  location  is 
tenned  ‘quarter’).  An  exam  pie  of  this  sch  erne  of  specific  elem  ents  used  to  calcu  late 
stress/strain  is  shown  in  Figure  3  for  a  10  by  10  element  mesh,  although  actual  composite 
plate  mesh  is  finer. 

The  normalized  transient  responses  of  di  splacement  and  strain  energy  typically 
show  the  same  shape  and  frequency,  with  onl  y  minor  differences  in  relative  am  plitudes, 
and  thus  can  be  used  interchangeably  to  demonstrate  the  behavior  of  the  composite  plate. 
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Figure  3.  Sample  of  Element  Locations  used  to  Calculate  Stress/Strain 

B.  DYNAMIC  RESPONSE  OF  COMPOSITE  PLATE  SUBJECTED  TO 
CONCENTRATED  FORCE  AND  CLAMPED  BOUNDARY 

The  baseline  stiffened  composite  plate  of  density  2020  kg/m3  and  elastic  modulus 
of  1 ,7e  10  P  a  with  clamped  edges  and  centrally  applied  concen  trated  force  w  ill  be 
discussed  first.  Folio  w  on  sectio  ns  will  exam  ine  variations  in  boundary  condition  , 
loading,  size,  shape  and  im  pact.  Figures  4,  5  and  6  show  the  response  of  the 
displacement,  strain  energy  and  kinetic  energy  of  the  plate  respectivel  y  (the  dry  case  is 
used  for  normalization). 


Figure  4. 


Normalized  Displacement  at  Center  of  Top  Skin  Plate 
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Figure  5.  Normalized  Strain  Energy  of  Composite  Plate 


Figure  6.  Normalized  Kinetic  Energy  of  Composite  Plate 

These  figures  show  the  com  parison  be  tween  one-side  wet  and  two-sides  wet 
structural  responses.  The  FSI  with  either  one-side  or  tw  o-sides  wet  of  the  com  posite 
structure  significantly  influences  both  the  m  agnitude  and  frequency  of  the  strain  energy 
plot.  The  oscillating  m  agnitude  and  the  fir  equency  are  drastically  reduced  by  the  FSI 
effect.  Two-sides  wet  FSI  results  in  the  lowe  st  peak  energy  values  and  their  frequency 
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among  the  three  cases.  However,  the  m  agnitude  of  oscillatory  behavior  is  the  least  for 
the  one-s  ide  wet  stru  cture.  The  figures  sh  ow  the  effects  of  FSI,  with  average 
displacement  and  energy  being  reduced  thr  ough  the  fluid  in  teraction.  Ad  ditionally,  FSI 
causes  a  decrease  in  frequency  and  m  agnitude  of  structural  responses,  with  significantly 
more  rapid  damping  effects  than  the  dry  case. 

The  transverse  displacem  ent  plot  at  the  node  of  the  applied  force  is  compared  in 
Figure  4  for  three  different  cases.  The  displace  ment  response  is  very  sim  ilar  to  that  of 
strain  energy  of  Figure  5.  The  two-sides  wet  structure  has  the  lowest  peak  displacem  ent 
and  frequency,  and  the  one-side  w  et  structur  e  has  the  least  vibrat  ory  m  otion.  It  is 
interesting  to  note  that  even  though  the  disp  lacement  characteristics  are  quite  different 
among  the  three  conditions;  their  respective  average  values  are  comparable. 

When  average  values  of  the  th  ree  strain  energy  variations  are  com  pared  (Figure 
5),  the  dry  structure  has  the  greatest  average  value  and  the  two-sides  wet  structure  has  the 
smallest  value.  Furthermore,  the  two-side  s  wet  structure  shows  energy  dissipation  as  a 
function  of  time. 

As  the  kinetic  energy  of  the  stiffened  st  ructure  is  compared  under  three  different 
surrounding  m  edia,  as  shown  in  F  igure  6,  the  dry  structure  show  s  a  very  significant 
oscillatory  behavior.  O  n  the  other  hand,  the  oscillation  of  kinetic  energy  is  suppressed 
quickly  for  the  wet  cases.  The  kinetic  energy  of  the  two-sides  wet  structure  is  the  lowest. 
The  two-sides  wet  structure  displays  the  fastest  decay  rate  of  the  kinetic  energy. 

The  normal  and  shear  strains  for  th  is  clamped  case,  for  each  of  the  locations  of 
interest  (center,  side  and  quarter),  are  shown  in  Figure  7,  with  normalization  with  respect 
to  dry  plate  x-axis  normal  strain.  Com  parison  of  the  nonn  al  strain  along  the  x-axis  also 
indicates  reduced  strain  s  for  wet  structures.  Wet  structures  have  very  high  frequency 
components  in  the  strain  response.  Howeve  r,  the  base  frequencie  s  of  both-side  wet 
structures  are  clearly  shown  lower  than  those  of  the  dry  structure.  Average  strain  v  alues 
are  more  or  less  sim  ilar  even  though  the  dry  structure  has  greater  am  plitudes  of  strain 
oscillation. 
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Norm. E>S, ,  Qrtr  Bottom  Plate  Norm. E>!x,  Side  Bottom  Plate  Norm.e  ,  Cntr  Bottom  Plate 


Figure  7. 


Normal  and  Shear  Strains  for  Clamped  Boundary  with  Concentrated  Force 
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C.  CLAMP  VERSUS  SIMPLE  SUPPO  RT  B  OUNDARY  CONDITION  WITH 

CONCENTRATED  FORCE 

Comparison  of  cla  mped  versus  simple  boundary  shows  little  difference  in 
dynamic  response  for  concentrated  force  load  ing.  The  displacem  ent  and  strain  energy 
plots  are  shown  in  Appendix  A.  The  kinetic  energy  responses  shown  in  Figures  8,  9  and 
10  shows  the  com  parison  for  the  dry,  two-sides  wet  and  one-side  wet  structures 
respectively.  There  is  aim  ost  no  differen  ce  for  t  he  dr  y  st  ructure;  however  ,  t  he  wet 
structures  show  slight  increase  in  energy.  This  is  expected  due  to  the  increased  degree  of 
freedom,  although  the  increased  energy  is  not  significant. 


Figure  8.  Comparison  of  Kinetic  Energy  of  Dry  Structure  between  Clamped  and 

Simple  Boundary 


Figure  9.  Comparison  of  Kinetic  Energy  of  Two-sides  Wet  Structure  between 

Clamped  and  Simple  Boundary 
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Figure  10.  Comparison  of  Kinetic  Energy  of  One-side  Wet  Structure  between 

Clamped  and  Simple  Boundary 

Similarly,  the  strains  at  the  center  an  d  side  locations  are  nearly  identical  and  are 
shown  in  Appendix  A.  Of  interest  are  the  quarter  location  strains,  which  show  som  e 
variance  between  the  boundary  condition  types,  with  the  clamp  condition  having  slightly 
higher  strains  for  the  dry  and  wet  structures  as  shown  in  Figure  11.  The  increase  in  strain 
for  the  clamped  boundary  was  expected  due  to  restricted  degree  of  freedom;  however  it  is 
surprising  to  be  evident  at  only  the  quarter  location. 

With  an  applied  conce  ntrated  f  orce,  th  ere  is  little  d  ifference  betwe  en  the  two 
types  of  boundary  conditions,  clamped  or  simple  support.  While  there  is  minor  increase 
in  kinetic  energy  of  the  wet  cas  es  for  simple  support  and  m  inor  decrease  in  strain  at  the 
quarter  location  of  the  com  posite  plate  for  the  sim  pie  support,  it  is  not  significant.  The 
FSI  effects  are  consistent  between  the  two  boundary  conditions. 
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Norm,  g  ,  Qrtr  Bottom  Plate,  1  Wet  Norm,  e  ,  Qrtr  Bottom  Plate,  2Wet  Norm,  e  ,  Qrtr  Bottom  Plate,  Dry 

XX  XX  XX  3 


Figure  1 1 .  Normal  and  Shear  Strain  Comparison  at  Quarter  Position  for  Clamped 
versus  Simple  Boundary  with  Concentrated  Force 
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D. 


CONCENTRATED  F  ORCE  VERSUS  P  RESSURE  LOADING  WITH 
CLAMPED  BOUNDARY 


Next,  the  dynamic  response  of  thin  composite  plate  was  compared  under  different 
loading  con  ditions:  con  stant  concentrated  fo  rce  and  equiv  alent  uniform  pressure,  each 
with  clamped  boundary.  The  basis  for  com  parison  is  clam  ped  boundary  with  constant 
concentrated  force  of  1000  N  applied  at  center  of  plate.  The  equivalent  unifonn  pressure 
loading  is  detenn  ined  from  the  concentrated  force  being  uniform  ly  a  pplied  over  the 
surface  of  the  0.3m  by  0.3m  plate,  giving  a  unifonn  pressure  load  of  1 1,1 1 1  Pa. 

The  comparison  for  the  dry  structure  under  the  two  loading  conditions  is  shown  in 
Appendix  B.  Under  dry  conditions,  the  pre  ssure  loading  versus  concentrated  force 
increases  the  amplitude  of  oscillation  for  displacement,  strain  energy,  and  kinetic  energy 
with  no  shift  in  frequency.  The  strain  at  the  center  location  has  increased  amplitude  but 
lower  average  strain.  T  he  normal  average  stra  in  at  the  side  location  is  incre  ased,  while 
the  shear  strain  is  com  parable  between  the  two  loading  conditions.  The  quarter  location 
exhibits  similar  strain  behavior  for  applied  force  and  pressure  loading. 

The  dry  structure  is  used  to  nonn  alize  the  wet  structure  responses  and  the 
displacement  response  showing  the  FSI  effects  are  shown  in  Figure  12.  The  wet  structure 
comparison  of  strain  and  kinetic  energy  for  force  versus  pressure  load  is  shown  in  Figures 
13  and  14,  respectively. 


Figure  12.  Wet  Structure  Displacement  Comparison  between  Force  and  Pressure 

Loading  with  Clamped  Boundary 


17 


1 


Figure  13.  Wet  Structure  Strain  Energy  Comparison  between  Force  and 
Pressure  Loading  with  Clamped  Boundary 


Figure  14.  Wet  Structure  Kinetic  Energy  Comparison  between  Force  and 
Pressure  Loading  with  Clamped  Boundary 

Figures  12  to  14  show  t  he  comparison  between  one-side  wet  and  two-sides  wet 

structural  responses.  The  FSI  with  both  one-side  and  two-sides  wet  reduces  the 

oscillating  magnitude  and  frequency  of  the  response  over  dry  structure.  The  pressure 

load  tends  to  produce  larger  am  plitude  of  os  cillation  than  concentrated  force,  but  th  e 

average  energy  is  similar,  while  the  mean  displacement  under  pressure  load  is  less.  Two- 

sides  wet  FSI  results  in  the  lowest  peak  energy,  peak  displacement,  and  frequency  among 

the  three  cases.  However,  the  m  agnitude  of  oscillatory  behavior  is  the  least  for  the  two- 

sides  wet  structure.  The  figures  sho  w  the  effects  of  FSI,  wi  th  average  displacement  and 
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energy  being  reduced  through  the  fluid  intera  ction.  Additionally,  FSI  causes  a  decrease 
in  frequency  and  m  agnitude  of  structural  responses  with  significantly  m  ore  rapid 
damping  effects  than  the  dry  case.  Of  note  in  Figure  14  is  the  slower  initial  response  of 
kinetic  energy  under  pressure  load. 

The  strain  behavior  for  the  dry  structure  is  shown  in  Appendix  B.  The  pressure 
load  vice  concentrated  force  comparison  show  increased  amplitude  of  strain  oscillation  at 
center  location  and  reduced  average  strain  as  well.  The  strain  response  at  the  side  and 
quarter  locations  was  similar  for  both  normal  and  shear  strains,  with  exception  of  norm  al 
strain  at  the  side  location  ha  ving  a  slightly  higher  m  agnitude  under  applied  pressure  than 
applied  force.  The  strain  behavior  f  or  wet  structures,  also  shown  in  Appendix  B,  was 
similar  to  that  of  dry,  with  an  increased  amplitude  oscillation  at  center  location,  but  with 
reduced  average  strain.  The  side  and  quarter  location  strain  response  of  wet  structures 
also  follow  ed  that  of  dry,  with  sam  e  diff  erences  of  the  norm  al  side  location  s  train 
exhibiting  higher  m  agnitude  under  pressure  lo  ading.  This  m  eans  the  concentrated  force 
has  a  greater  FSI  effect  than  the  pressure  at  the  side  and  quarter  locations. 

E.  CONCENTRATED  F  ORCE  VERSUS  P  RESSURE  LOADING  WITH 

SIMPLE  SUPPORT  BOUNDARY 

Next  the  dynam  ic  response  of  a  thin  com  posite  plate  was  com  pared  under 
different  loading  conditions,  constant  con  centrated  force  and  equivalent  unif  orm 
pressure,  with  a  sim  pie  support  boundary.  The  basis  for  com  parison  is  sim  pie  support 
boundary  with  constant  concentrated  force. 

The  com  parison  for  the  dry  structu  re  under  th  e  two  loading  conditio  ns  with 
simple  support  is  shown  in  Appendix  C.  Under  dry  conditions  the  pressure  and 
concentrated  force  loading  have  nearly  identi  cal  responses  with  no  di  scemable  change  in 
amplitude  or  frequency  for  displacement,  strain  energy  and  kinetic  energy.  The  strains  at 
the  center  locations  have  nearly  identical  response  for  applied  force  and  pressure  loading. 
The  norm  al  averag  e  strains  at  th  e  side  and  qu  arter  lo  cation  is  in  creased  for  press  ure 
loading,  wh  ile  the  shear  strain  sho  ws  higher  amplitude  of  oscillation  for  the  pres  sure 
loading  condition. 
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The  wet  structure  response  com  parison  with  simple  boundary  for  force  versus 
pressure  load  is  shown  in  Figures  15,  16  and  17  for  displacem  ent,  strain  energy  and 
kinetic  energy  respectively. 


Figure  15.  Wet  Structure  Displacement  Comparison  of  Force  and  Pressure  Loading 

with  Simple  Support  Boundary 


Figure  16.  Wet  Structure  Strain  Energy  Comparison  of  Force  and  Pressure  Loading 

with  Simple  Support  Boundary 
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Figure  17.  Wet  Structure  Kinetic  Energy  Comparison  of  Force  and  Pressure  Loading 

with  Simple  Support  Boundary 

The  FSI  with  both  one-side  and  two-sides  wet  structures  reduce  the  os  cillating 
magnitude  and  frequency  of  the  response  over  the  dry  structure.  With  a  simple  boundary, 
the  pressure  and  force  load  track  very  well  with  one  another,  with  only  minor  difference 
in  frequency  evident  in  the  two-sides  wet  structure  displacem  ent  and  strain  energy. 
Unlike  the  clam  ped  boundary,  there  is  no  dela  y  in  response  of  kinetic  energy  with  a 
simple  boundary  for  the  force  and  pressure  loa  d.  Again,  the  two-side  s  wet  FSI  results  in 
the  lowest  peak  energy,  peak  displacem  ent  and  frequency  among  th  e  three  cases.  The 
figures  sho  w  the  effects  of  FSI,  with  aver  age  displacem  ent  and  energ  y  being  red  uced 
through  the  fluid  interaction. 

The  strain  behaviors  for  the  three  stru  ctures  are  shown  in  Appendix  C.  The 
pressure  and  force  load  strains  track  each  other  using  a  sim  pie  support  boundary  at  the 
center  location.  The  strain  respon  se  at  th  e  side  location  has  sim  ilar  am  plitude  of 
oscillation,  with  the  norm  al  st  rains  slightly  higher  under  pr  essure  load.  The  strain 
behavior  at  the  quarter  location  is  sim  ilar  for  wet  structures,  although  the  two-sides  wet 
structure  has  less  amplitude,  the  wet  structures  overall  have  approximately  equal  average 
strain. 

F.  SIZE  OF  COMPOSITE  PLATE 

Next  the  influence  of  com  posite  plate  size  on  FSI  is  exam  ined  by  increas  ing  the 

size  of  the  square  plate  from  0.3m  to  0.5m  on  a  side.  The  com  parison  is  m  ade  using 
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clamped  boundary  condition  with  applied  concentrated  force.  The  displacement  response 
for  the  three  structures  is  shown  in  Appendix  D,  and  indicates  that  increases  in  plate  size 
yield  a  decrease  in  frequency,  with  the  tw  o-sides  wet  structure  having  a  substantial 
decrease  in  frequency.  Also,  FSI  damping  is  slower  as  the  plate  size  in  creases.  Similar 
results  are  visible  in  strain  and  kinetic  en  ergy  response  betw  een  the  two  sizes  of  plates 
shown  in  Figures  18,  19  and  20,  for  the  dry,  two-sides  wet  and  one-s  ide  wet  structures 
respectively. 


Figure  18.  Dry  Structure  Strain  and  Kinetic  Energy  Comparison  for  Plate  Size 
Variations  with  Concentrated  Force  and  Clamped  Boundary 


Figure  19.  Two-sides  Wet  Structure  Strain  and  Kinetic  Energy  Comparison  for  Plate 
Size  Variations  with  Concentrated  Force  and  Clamped  Boundary 
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Figure  20.  One-side  Wet  Structure  Strain  and  Kinetic  Energy  Comparison  for  Plate 
Size  Variations  with  Concentrated  Force  and  Clamped  Boundary 

The  larger  size  plate  has  lower  frequenc  y,  slower  initial  response,  and  slower 
long-term  damping.  The  FSI  of  the  larger  pi  ate  is  les  s  and  thus  has  higher  am  plitude 
oscillations.  The  dry  structure  has  similar  average  energies  between  the  two  plate  sizes, 
while  the  steady  state  energy  of  the  larger  plate  is  marginally  greater  for  the  one-side  wet 
structure.  The  difference  in  energy  between  the  two  plate  sizes  is  more  pronounced  in  the 
two-sides  wet  structure,  where  the  kinetic  energy  clearly  shows  th  e  significant  delay  in 
response  due  to  the  damping  effect  of  fluid. 

The  comparison  of  strain  between  the  two  sized  plates  is  shown  in  Appendix  D. 
However,  there  is  no  clearly  identifiable  characteristic  between  the  strains  with  exception 
of  some  decreased  frequency  and  comparable  average  normal  and  shear  strains. 

G.  SHAPE  OF  COMPOSITE  PLATE 

The  influence  of  composite  plate  shape  on  FSI  is  exam  ined  next  by  eh  anging  the 

shape  of  the  plate  from  square  to  rectangular  while  m  aintaining  equivalent  area,  thus  the 

rectangular  plate  is  0.2  m  by  0.45m.  As  with  com  parison  of  plate  size,  the  shape 

comparison  is  made  using  clamped  boundary  conditions  with  applied  concentrated  force. 

The  displacement  response  for  the  three  structures  is  shown  in  Appendix  E,  and  indicates 

the  rectangular  shape  has  increase  in  frequency  and  decrease  in  amplitude  of  oscillation 

over  the  square  plate  of  equivalent  area,  with  the  average  displacem  ent  of  the  th  ree 

structures  (dry,  two-sides  wet,  one-side  wet)  slightly  greater  for  the  rectangular  shape. 
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The  strain  and  kine  tic  energy  response  between  the  two  shap  es  of  plates  are 
shown  in  Figures  21,  22  and  23,  for  the  dry,  two-sides  wet  and  one-s  ide  wet  structures 
respectively.  The  rectangular  plate  has  a  higher  frequency  and  faster  damping  rate.  The 
rectangular  plate  has  lo  wer  a  mplitude  of  oscillations.  The  average  energies  of  the 
rectangular  plate  are  higher  th  an  those  of  the  square  plat  e.  The  difference  in  energy 
between  the  two  plate  shapes  is  m  ore  pronounc  ed  in  the  two-sides  wet  structure,  w  hich 
clearly  shows  the  FSI  effect  is  greates  t  for  two-sides  wet  structur  e  and  the  overall  FSI 
effect  is  less  for  the  rectangular  vice  square  plate  since  the  peak  energy  is  greater. 


Figure  2 1 .  Dry  Structure  Strain  and  Kinetic  Energy  Comparison  for  Plate  Shape 
Variations  with  Concentrated  Force  and  Clamped  Boundary 


Figure  22.  Two-sides  Wet  Structure  Strain  and  Kinetic  Energy  Comparison  for  Plate 
Shape  Variations  with  Concentrated  Force  and  Clamped  Boundary 
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Figure  23.  One-side  Wet  Structure  Strain  and  Kinetic  Energy  Comparison  for  Plate 
Shape  Variations  with  Concentrated  Force  and  Clamped  Boundary 

The  comparison  of  strain  between  the  two  shaped  plates  is  shown  in  Appendix  E. 
The  shear  strains  for  the  rectangular  plate  are  all  less  than  for  the  e  quivalent  area  square 
plate.  The  average  norm  al  strain  at  center  lo  cation  of  all  three  structures  is  a  little  more 
for  the  rectangular  vice  the  squa  re  plate,  while  the  nonn  al  strains  at  the  side  and  quarter 
locations  are  very  sim  ilar.  Overall,  the  re  duction  in  energy  due  to  FSI  effects  of  the 
rectangular  plate  shape  is  less  than  the  square  plate. 

H.  COMPOSITE  DENSITY 

Next  the  influence  of  com  posite  m  aterial  density  on  dynam  ic  response  is 
examined.  Since  the  response  of  displacement  is  similar  to  strain  energy,  only  the  strain 
energy  will  be  used  here  and  displacem  ent  plots  are  in  Appendix  F.  Figures  24,  25  and 
26  show  the  strain  and  ki  netic  energy  for  the  dry,  two-  sides  wet  and  one-side  wet 

■3 

structures  respectively,  w  ith  each  u  sing  a  composite  plate  of  density  2020  kg/m  and 
elastic  m  odulus  of  1 ,7e  10  Pa,  with  concentrated  force  and  clam  ped  boundary  for 
normalization.  For  the  dry  structu  re  it  is  clearly  visible  that  increasing  density  causes  a 
decrease  in  frequency,  however,  due  to  FSI  this  feature  is  not  as  pronounced  in  the  w  et 
structures.  The  wet  structures  show  only  s  light  difference  in  frequency  and  the  peak 
strain  energ  y  occurs  in  lowest  density  with  on  ly  m  inimal  decrease  in  peak  energy  as 
density  increases.  The  kineti  c  energy  shows  a  faster  rate  of  da  mping  with  increasing 
density  for  the  wet  structures. 
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Figure  24.  Dry  Structure  Strain  and  Kinetic  Energy  Comparison  for  Density 
Variations  with  Concentrated  Force  and  Clamped  Boundary 


Figure  25.  Two-sides  Wet  Structure  Strain  and  Kinetic  Energy  Comparison  for 
Density  Variations  with  Concentrated  Force  and  Clamped  Boundary 
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Figure  26.  One-side  Wet  Structure  Strain  and  Kinetic  Energy  Comparison  for 
Density  Variations  with  Concentrated  Force  and  Clamped  Boundary 


To  highlight  the  specific  FSI  effects,  the  m  ethod  of  normalization  was  altered 
such  that  each  wet  structure  is  norm  alized  to  its  respective  dry  stru  cture  and  the  th  ree 
different  densities  are  plotted  together  to  show  its  effect  on  response.  W  ith  this  alternate 
nonnalization,  the  strain  and  ki  netic  energy  is  shown  in  Fi  gures  27  and  28  for  the  two- 
sides  wet  and  one-side  wet  struct  ures,  respectively.  In  this  representation,  it  is  clear  FSI 
gives  a  reduction  in  peak  energy,  is  m  ore  signif  icant  in  two-sid  es  wet  structu  re  and 
drastically  reduces  the  high  frequency  oscillat  ion  from  the  dry  structure.  Also,  the  two- 
sides  wet  structure  kinetic  energy  shows  faster  response,  as  density  increases,  from  initial 
load  application  to  peak  energy  value  and  subsequent  decay  toward  steady  state  value. 


Figure  27.  Two-sides  Wet  Structure  Strain  and  Kinetic  Energy  Comparison  for 
Density  Variations  with  Concentrated  Force  and  Clamped  Boundary  (Alternate 

Normalization) 
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Figure  28.  One-side  Wet  Structure  Strain  and  Kinetic  Energy  Comparison  for 
Density  Variations  with  Concentrated  Force  and  Clamped  Boundary  (Alternate 

Normalization) 


The  variations  in  strain  of  the  three  structures  are  shown  in  Appendix  F.  For  the 
dry  and  one-side  wet  structure,  a  decreas  e  in  frequency  as  d  ensity  increases  is  evid  ent, 
but  is  not  identif  iable  f  or  the  two-  sides  we  t  structure.  N  one  of  the  structures  exhibit 
significant  variation  in  magnitude  of  strain  for  the  different  density  values.  As  shown  i  n 
Figure  29,  using  the  previously  discussed  alternate  normalization,  the  strains  for  the  two- 
sides  wet  structure  have  sim  ilar  relative  magnitude,  for  the  three  loca  tions  (center,  side, 
quarter),  with  no  discemable  shift  in  frequency  due  to  density  variations,  while  Figure  30 
shows  the  one-side  wet  structure  having  a  deer  ease  in  f  requency  from  density  inc  reases 
with  similar  relative  strain  magnitude. 
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Norm.  E>s, ,  Qrtr  Bottom  Plate  Norm.  E)(x ,  Side  Bottom  Plate  Norm.  E  i  Cntr  Bottom  Plate 


Figure  29.  Normal  and  Shear  Strains  for  Comparison  of  Differ  Density  for  Two-sides 
Wet  Structure  with  Concentrated  Force  and  Clamped  Boundary  (Alternate 

Normalization) 
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Norm.Exx,Qrtr  Bottom  Plate  Norm. Ejk,  Side  Bottom  Plate  Norm.e  ,  Cntr  Bottom  Plate 


Figure  30.  Normal  and  Shear  Strains  for  Comparison  of  Differ  Density  for  One-side 
Wet  Structure  with  Concentrated  Force  and  Clamped  Boundary  (Alternate 

Normalization) 
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I.  COMPOSITE 


MODULUS 


The  influence  of  com  posite  m  aterial  elastic  modulus  on  dynam  ic  response  is 
examined  next.  The  d  isplacement  response  fo  r  different  elastic  m  odulus  is  shown  in 
Appendix  G.  Figures  3 1,  32  and  33  show  the  strain  and  kinetic  energy  for  the  dry,  two- 
sides  wet  and  one-side  we  t  structures  respectively,  with  each  using  a  com  posite  plate  of 
density  2020  kg/m  and  elastic  m  odulus  of  1.7e  Pa,  with  concentrated  fore  e  and 
clamped  boundary  for  nonn  alization.  For  the  thr  ee  structures,  an  increase  in  frequency 
and  decrease  in  am  plitude  is  clearly  visi  ble  f  or  increasing  elastic  m  odulus.  As  the 
composite  elastic  m  odulus  increases,  the  stru  cture  becomes  stiffer  and  as  the  strain  and 
kinetic  energy  plots  show,  the  average  energy  decreases  with  increasing  modulus.  Also, 
the  amplitude  of  oscillation  decreases  with  increasing  modulus.  The  wet  structures  show 
a  similar  rate  of  damping  with  increasing  modulus. 


Figure  3 1 .  Dry  Structure  Strain  and  Kinetic  Energy  Comparison  for  Elastic  Modulus 
Variations  with  Concentrated  Force  and  Clamped  Boundary 
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Figure  32.  Two-sides  Wet  Structure  Strain  and  Kinetic  Energy  Comparison  for 
Elastic  Modulus  Variations  with  Concentrated  Force  and  Clamped  Boundary 
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Figure  33.  One-side  Wet  Structure  Strain  and  Kinetic  Energy  Comparison  for  Elastic 
Modulus  Variations  with  Concentrated  Force  and  Clamped  Boundary 

To  highlight  the  specific  FSI  effects,  the  m  ethod  of  normalization  was  altered  as 
discussed  in  the  previous  section.  With  this  alternate  normalization,  the  strain  and  kinetic 
energy  is  shown  in  Figures  34  and  35  for  the  two-sides  wet  and  one-side  wet  structures, 
respectively.  In  this  representation,  it  is  cl  ear  F  SI  gives  a  reduction  in  peak  energy,  is 
more  significant  in  two-sides  wet  s  tructure  and  drastically  reduces  the  high  frequency 
oscillation.  The  average  strain  energy  is  co  mparable  at  each  of  the  different  elastic 
modulus  values,  with  only  the  frequency  and  amplitude  varying,  while  the  kinetic  energy 
tends  to  decrease  with  increasing  modulus. 
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Figure  34.  Two-sides  Wet  Structure  Strain  and  Kinetic  Energy  Comparison  for 
Elastic  Modulus  Variations  with  Concentrated  Force  and  Clamped  Boundary 

(Alternate  Normalization) 


Figure  35.  One-side  Wet  Structure  Strain  and  Kinetic  Energy  Comparison  for  Elastic 
Modulus  Variations  with  Concentrated  Force  and  Clamped  Boundary  (Alternate 

Normalization) 

The  strain  variations  of  the  three  structures  are  shown  in  Appendix  G.  For  these 
structures,  the  decrease  in  amplitude  of  oscillation  as  elastic  modulus  increases  is  evident 
while  the  in  crease  in  f  requency  les  s  noticea  ble  in  the  strain  plots.  The  m  ost  notable 
feature  is  the  large  re  duction  in  strain  as  the  m  odulus  increases;  this  is  due  to  the 
increased  stiffness.  Us  ing  the  alternate  norm  alization,  so  that  spec  ific  influences  of 
elastic  modulus  and  FS  I  can  be  highlighted  fo  r  strain,  the  strain  responses  of  the  two- 
sides  wet  and  one-side  wet  structures  are  shown  in  Figures  36  and  37  ,  respectively.  For 
the  two-sides  wet  structure,  th  e  relative  magnitude  of  strain  is  consistent  for  each  of  the 
three  moduli  across  the  three  different  locations  on  the  plate,  with  minor  indication  of  the 
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increase  in  frequency  with  increasing  m  odulus.  The  frequency  shift  is  m  uch  more 
evident  in  the  one  -side  wet  s  tructure,  whil  e  again  the  rela  tive  stra  in  m  agnitude  is 
consistent  for  the  various  modulus  values. 
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Time  (sec.) 


Time  (sec.) 


Time  (sec.) 


Figure  36.  Normal  and  Shear  Strains  for  Comparison  of  Differ  Elastic  Modulus  for 
Two-sides  Wet  Structure  with  Concentrated  Force  and  Clamped  Boundary 

(Alternate  Normalization) 
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Figure  37.  Normal  and  Shear  Strains  for  Comparison  of  Differ  Elastic  Modulus  for 
One-side  Wet  Structure  with  Concentrated  Force  and  Clamped  Boundary 

(Alternate  Normalization) 
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J.  IMPACT 


LOADING 


The  final  dynam  ic  behavior  exam  ined  is  the  impact  response  of  com  posite  plate 
from  a  steel  projectile.  Three  velocities  are  examined  for  the  projectiles:  1  m/s,  5  m/s  and 
10  m/s.  In  addition,  the  response  due  to  tw  o  im pact  face  shapes  are  com  pared,  each 
having  the  same  surface  area  for  i  mpact  and  equal  mass.  A  cylindrical  shaped  impacto  r 
has  a  circular  shape  area  of  impact  and  a  rectangular  shaped  impactor  has  a  square  shape 
area  of  impact.  Each  projectile  contacts  the  composite  plate  at  the  center. 

1.  Shape  of  Impactor 

To  investigate  any  dependence  on  shape  of  the  impact  object,  a  fixed  velocity  of 
10  m/s  is  used  to  com  pare  the  difference  in  response  between  the  tw  o  shapes  of  im  pact 
projectiles.  Figures  38,  39  and  40  compare  the  displacement,  strain  energy  and  kinetic 
energy  response  of  the  three  structu  res  to  the  circular  and  square  shape  impactor.  The 
figures  show  FSI  gives  a  significant  reduction  in  amplitude  and  frequency  and  the  square 
impact  face  has  less  am  plitude  for  dry  and  one-side  wet  s  tructure  than  the  circula  r  face 
impactor.  The  two-sides  wet  structure  show  s  sim  ilar  initial  response  between  the  two 
shapes  of  impact  and  then  the  square  face  impactor  respo  nse  stead  ies  out  with  higher 
amplitude  o  f  oscillation  .  The  average  di  splacement  is  co  mparable  b  etween  th  e  two 
impact  shapes,  while  the  average  energy  is  less  for  the  square  face  impact  for  each  of  the 
three  structures. 


Circular  Impact  Area 


Square  Impact  Area 


Figure  38.  Displacement  Comparison  of  Three  Structures  Due  to  Different  Impactor 

Shapes  at  10  m/s 
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Figure  39.  Strain  Energy  Comparison  of  Three  Structures  Due  to  Different  Impactor 

Shapes  at  10  m/s 
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Figure  40.  Kinetic  Energy  Comparison  of  Three  Structures  Due  to  Different  Impactor 

Shapes  at  10  m/s 


The  norm  al  strain  at  each  lo  cation  is  com  pared  in  Figure  41  for  the  two 
shapes  of  impactor.  The  FSI  decreas  es  si  ightly  th  e  str  ain  am  plitude  with  m  inor 
decrease  in  frequency.  The  two-sides  wet  st  ructure  has  a  larger  peak  strain  for  the 
center  and  side  locations  for  th  e  square  faced  itn  pactor  over  the  circular  faced 
impactor,  while  the  cylindrical  impactor  ha  s  peak  strain  in  dry  and  one-side  wet 
structure.  The  shear  strain  is  slightly  higher  for  square  face  im  pactor  at  center 
location,  while  the  shear  strains  are  com  parable  between  im  pactor  at  the  side  and 
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quarter  locations.  The  FSI  is  m  ore  pronounced  for  the  square  impactor  one-side 
wet,  and  more  for  the  cylindrical  impactor  two-sides  wet  case. 
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Figure  4 1 .  Normal  Strain  Comparison  of  Three  Structures  Due  to  Different  Impactor 

Shapes  at  10  m/s 
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Circular  Impact  Area 


Square  Impact  Area 


Figure  42.  Shear  Strain  Comparison  of  Three  Structures  Due  to  Different  Impactor 

Shapes  at  10  m/s 
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Additional  figures  comparing  the  response  between  the  square  and  circular  faced 
impactor  are  contained  in  Appendix  H.  These  additional  figures  highlight  the  FSI 
differences  for  the  two  different  shapes  of  impactor  for  the  two-sides  wet  and  one-side 
wet  structures  by  norm  alizing  each  to  their  respective  dr  y  structure  response.  To 
summarize,  FSI  slightly  decrea  ses  the  strain  amplitude  and  frequency  of  oscillation  with 
the  square  impactor  having  a  slight  increase  in  frequency  and  am  plitude  of  oscillation 
over  the  circular  impactor  for  both  two-sides  and  one-side  wet.  The  square  impactor  also 
has  slightly  less  av  erage  energy  (strain  and  ki  netic).  The  stra  ins  a  re  nearly  the  s  ame 
between  the  cylinder  and  square  im  pactor  with  the  square  having  higher  peak  strain  at 
center  position  and  comparable  for  the  side  and  quarter  positions.  The  average  strains  are 
roughly  the  same  except  for  the  center  pos  ition  which  is  higher  due  to  higher  peak  strain 
initially. 

2.  Velocity  of  Impact 

The  effect  of  impact  velocity  is  straight  forward;  increasing  impact  velocity  gives 
increased  magnitude  of  plate  displacement,  strain  and  kinetic  energies.  When  combining 
the  varying  im  pact  velocities  with  different  s  haped  im  pactors,  there  are  som  e  slight 
differences  in  response.  The  shift  in  initial  response  when  comparing  the  three  velocities 
is  due  to  th  e  time  difference  required  for  the  impactor  to  trave  rse  the  distance  to  the 
composite  plate  and  should  not  be  misinterpreted  as  a  frequency  shift. 

The  response  of  each  o  f  three  structures  to  different  initial  im  pact  velocities  is 
shown  in  A  ppendix  I  for  both  circular  and  square  faced  impactors.  Increas  ing  impact 
velocity  s  imply  increas  es  the  respon  se.  Ge  nerally,  the  square  faced  impactor  has  less 
amplitude  of  oscillation  and  average  values  for  the  dry  and  one-side  wet  structure  and  the 
two-sides  wet  structure  amplitude  of  oscillation  and  average  value  is  s  imilar  for  the  two 
different  impactors.  The  normal  strains  are  comparable  with  only  very  slight  decrease  in 
amplitude  o  f  oscillation  for  the  square  impactor  for  each  o  f  the  th  ree  structures.  The 
shear  strain  is  also  similar  am  ong  the  three  structures  and  two  im  pactors  for  the  three 
impact  velocities. 
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To  focus  on  the  FSI  effects,  each  o  f  the  impact  velocities  for  the  two-s  ides  and 
one-side  wet  stru  ctures  are  norm  alized  to  a  respective  dry  structure.  T  hese  nonnalized 
responses  are  shown  in  Figures  43,  44  and  45  fo  r  displacement,  strain  energy  and  kinetic 
energy  respectively.  Using  this  normalization,  it  is  clear  FSI  causes  significant  decreases 
in  frequency  and  a  mplitude  range  of  res  ponse.  The  two-sides  wet  structure  shows 
decreased  peak  values  while  the  on  e-side  wet  has  slightly  increased  peak  values  over  the 
strictly  dry  structure  du  e  to  effects  of  the  wa  ter  layer  on  one  side  of  plate.  The  square 
faced  impactor  has  higher  relative  amplitude  of  oscillation  for  displacement,  but  slightly 
lower  energies  than  those  of  circular  faced  impactor. 
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Figure  43.  Comparison  of  Displacement  Response  Due  to  Impact  Velocity  Effects  for 

Circular  and  Square  Faced  Impactor 
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Figure  44.  Comparison  of  Strain  Energy  Response  Due  to  Impact  Velocity  Effects  for 

Circular  and  Square  Faced  Impactor 


The  normal  and  shear  strains  using  the  nor  malization  to  highlight  the  F  SI  effects 
are  shown  in  Figures  46  and  47,  respectively.  These  figures  show  the  relative  magnitudes 
of  strain  f  or  all  three  velocitie  s  ar  e  sim  ilar  with  the  ex  ception  of  the  two-sides  wet 
structure  with  square  face  impactor  has  slightly  increased  peak  strain. 
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Figure  45.  Comparison  of  Kinetic  Energy  Response  Due  to  Impact  Velocity  Effects 

for  Circular  and  Square  Faced  Impactor 


43 


Circular  Square 

Figure  46.  Comparison  of  Normal  Strain  Due  to  Impact  Velocity  Effects  for  Circular 

and  Square  Faced  Impactor 
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Figure  47.  Comparison  of  Shear  Strain  Due  to  Impact  Velocity  Effects  for  Circular 

and  Square  Faced  Impactor 
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V.  NUMERICAL  MODEL  COMPARISON  TO  EXPERIMENT 


The  principal  focus  of  this  work  is  nu  merical  study  of  various  param  eters  which 
affect  the  dynam  ic  be  havior  of  com  posite  plates.  A  separate  study  conducted 
experimentally  exam  ines  the  behavior  of  dry  and  wet  plates  s  ubjected  to  impact. 
Numerical  and  experim  ental  studies  each  have  their  respective  advantages  and 
disadvantages  and  are  used  to  com  plement  each  other.  In  particular,  the  v  arious 
parametric  studies  conducted  in  this  wo  rk  were  only  possible  utilizing  numerical 
modeling.  Experimental  testing  is  limited  to  measuring  forces  and  strains  through  gages. 
Preliminary  comparison  of  experimental  and  numerical  work  is  noted  here  to  detennine 
methods  for  improvement  to  follow  on  research. 

A.  EXPERIMENTAL  SETUP  FOR  IMPACT  LOADING 

Preliminary  experimental  behavior  study  of  thin  composite  plates  is  conducted  on 
12  inch  square  and  1/16  inch  thick  plate  clamped  in  the  frame  of  an  impact  testing  device. 
The  device  uses  a  weighted  sled  system  to  strike  a  cylindrical  impactor.  Multiple  impacts 
are  prevented  using  a  large  sp  ring  opposing  the  cylindrical  impactor  such  that  only  one 
impact  event  takes  place.  There  is  a  for  ce  measuring  gage  m  ounted  on  the  end  of  the 
cylindrical  impactor  to  m  easure  the  force  durin  g  contact  with  the  com  posite  plate.  The 
schematic  of  the  experimental  device  setup  is  shown  in  Figure  48.  The  underside  of  the 
composite  plate  is  instrumented  with  strain  gages,  bonded  to  the  pi  ate  with  epoxy,  in  the 
layout  shown  in  Figure  49.  Th  e  strain  gages  measure  approxi  mately  1  cm  square.  Gage 
2  is  in  the  center  location  for  comparison  to  the  numerical  model  a  nd  is  directly  below 
impact  site.  Gages  1  and  4  are  repr  esentative  of  a  side  location  sim  ilar  to  the  num  erical 
model  and  gages  3  an  d  5  are  sim  ilar  to  the  quarter  location.  The  data  acqu  isition 
software  measures  the  transient  force  and  strain  data  at  1000  Hz  sampling  rate. 

B.  NUMERI  CAL  MODEL 

The  1/1 6  th  inch  thick,  12  inch  by  12  inch  com  posite  plate  is  m  odeled  using  shell 

elements  with  a  m  esh  seed  of  60  nodes  per  side  .  The  m  esh  size  is  chosen  to  adequately 

approximate  the  im  pact  force  gage  area  and  strain  gage  size  reasonab  ly.  The  i  mpact 
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force  m  easured  experim  entally  is  converted  to  equivalent  pressure  and  applied  to 
elements  approximating  a  cylindrical  impactor  striking  the  plate.  In  the  numerical  model, 
stress  is  co  mputed  and  strains  are  calcu  lated  using  standard  stress-strain  transform  ation 
equations.  The  strain  over  the  area  of  the  numbered  experim  ental  strain  gages  is 
calculated  b  y  averag  ing  the  elem  ents  which  ap  proximate  the  size  of  the  strain  gage  to 
compare  with  the  experimentally  measured  strains. 


(a)  Schematic 


(b)  As-built 


Figure  48.  Impact  Device  Experimental  Setup 
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(a)  Strain  Gage  Layout  (b)  Actual  Composite  Plate 

Figure  49.  Experiment  Strain  Gage  Layout  on  Underside  of  Composite  Plate 
(Dimensions  in  parenthesis  are  given  in  inches) 


C.  COMPARISON  OF  EXPERIMENTAL  AND  NUMERICAL  RESULTS 

One  of  the  challenges  experim  ental  work  presents  when  m  easuring  strains  is 
getting  a  good  bond  between  the  strain  gages  an  d  the  composite  plate.  Another  is  having 
the  strain  g  ages  aligne  d  perf  ectly  with  the  direction  of  fibers  in  the  com  posite  and 
minimizing  the  area  of  the  gage  covering  th  e  matrix  which  f  orms  the  com  posite.  The 
experiment  is  conducted  in  a  one-side  wet  sc  enario  using  an  anechoic  tank  to  m  inimize 
water  disturbance  effects,  with  the  side  oppos  ite  of  i  mpact  on  the  composite  plate  kept 
dry  through  a  plexi-glass  box  bound  to  the  underside  of  the  com  posite  plate.  The 
experiment  is  also  run  in  a  com  pletely  dry  condition.  Both  dry  and  wet  cases  use  the 
same  impact  force  by  dropping  the  weighted  sled  from  full  height,  giving  the  steel  impact 
rod  roughly  a  5  m  /s  initial  ve  locity.  The  experim  ents  were  first  conducted  in  the  wet 
condition  and  then  dry.  Following  the  wet  experim  ents  it  was  identified  that  the  strain 
gage  labeled  Gage  1  had  broken  free  from  the  com  posite  plate  and  hence  was  not 
available  for  the  dry  experiment. 

The  comparison  of  normal  strain  for  gage  1  location,  between  the  experim  ent  and 
simplified  Finite  Element  model  is  shown  in  Figure  50.  A  s  shown,  there  is  not  a  good 
comparison  between  the  experim  ent  and  m  odel  results  for  the  one-side  wet  condition 
(note  that  the  strain  gage  fell  off  prior  to  dry  experiment  and  there  is  no  experimental  data 
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to  compare).  Because  the  strain  gage  fell  off  after  the  wet  experiment,  the  data  shown  is 
possibly  erroneous  due  to  the  strain  gage  disbonding  and  any  com  parison  at  the  gage  1 
location  is  suspect. 


Dry  Condition 


One-side  Wet  Condition 


Figure  50.  Comparison  of  Normal  Strain  at  Gage  1  Location  Between  Experiment 

and  FEM  in  Dry  and  One-side  Wet  Condition 


The  comparison  of  normal  strain  for  gage  2  location,  between  the  experim  ent  and 
simplified  FEM  m  odel  is  shown  in  Figure  51.  As  shown,  the  com  parisons  between  the 
experiment  and  model  results  are  quite  good  for  both  the  one-side  wet  and  dry  condition. 
This  good  agreem  ent  between  the  experim  ent  and  numerical  model  is  evidence  of  the 
feasibility  to  accurately  pred  ict  co  mposite  plate  response  us  ing  finite  elem  ent  models. 
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This  gives  more  flexibility  for  researchers  as  many  more  parameters  can  be  varied  with  a 
numerical  model.  The  fact  that  the  x-axis  model  strain  is  higher  than  the  experiment  and 
the  y-axis  is  lower  is  an  indication  there  m  ay  be  som  e  misalignment  of  t  he  strain  gage 
with  the  fiber  direction.  If  this  is  the  case,  som  e  improvement  can  be  obtained  through 
use  of  a  Mohr  Circle  transformation. 


Dry  Condition 


One-side  Wet  Condition 


Figure  5 1 .  Comparison  of  Normal  Strain  at  Gage  2  Location  Between  Experiment 

and  FEM  in  Dry  and  One-side  Wet  Condition 


The  comparison  of  nonnal  strain  for  the  gage  3  location,  between  the  experim  ent 
and  simplified  FEM  model  is  shown  in  F igure  52.  As  shown,  the  com  parisons  between 
the  experiment  and  model  results  are  quite  good  for  the  dry  condition  but  not  the  one-side 
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wet  condition.  Use  of  Mohr  Circle  transf  ormation  m  ay  i  mprove  the  dry  com  parison. 
What  is  encouraging  is  the  trend  between  e  xperiment  and  model  tr  acks.  Unfortunately, 
there  is  not  a  good  explanation  of  why  the  dr  y  condition  is  in  such  good  agreem  ent  but 
the  one-side  wet  condition  is  not. 


Dry  Condition 


One-side  Wet  Condition 


Figure  52.  Comparison  of  Normal  Strain  at  Gage  3  Location  Between  Experiment 

and  FEM  in  Dry  and  One-side  Wet  Condition 


The  comparison  of  normal  strain  for  gage  4  location,  between  the  experim  ent  and 
simplified  FEM  m  odel  is  shown  in  Figure  53.  As  shown,  there  is  good  agreem  ent 
between  the  experiment  and  model  results  for  the  dry  cond  ition  and  the  trends  agree  for 
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the  one-side  wet  condition,  although  the  m  agnitudes  are  off.  Again  a  Mohr  Circle 
transformation  could  improve  the  dry  and  wet  comparison. 


Dry  Condition 


One-side  Wet  Condition 


Figure  53.  Comparison  of  Normal  Strain  at  Gage  4  Location  Between  Experiment 

and  FEM  in  Dry  and  One-side  Wet  Condition 


The  comparison  of  normal  strain  for  gage  5  location,  between  the  experim  ent  and 
simplified  FEM  m  odel  is  shown  in  Figure  54.  As  shown,  the  trends  between  the 
experiment  and  m  odel  are  sim  ilar,  but  the  magnitudes  are  not,  and  application  of  Mohr 
Circle  will  not  im  prove  the  values  as  both  the  x-axis  and  y-  axis  normal  strains  are  over 
predicted  in  the  numerical  model. 
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Dry  Condition 


One-side  Wet  Condition 


Figure  54.  Comparison  of  Normal  Strain  at  Gage  5  Location  Between  Experiment 

and  FEM  in  Dry  and  One-side  Wet  Condition 


In  summary,  the  strain  in  the  vicinity  of  impact,  at  gage  2  location,  compares  very 
well  between  model  and  experiment.  Moving  away  from  impact  location  either  to  side  or 
quarter  lo  cation  r  esults  in  le  ss  agr  eement  be  tween  the  experim  ent  and  m  odel.  This 
indicates  that  proper  strain  gage  alignm  ent  with  fiber  direction  and  good  bonding  over 
fiber  vice  matr  ix  is  im  portant.  Other  things  to  consider  in  f  uture  work  are  altering  the 
element  size  in  the  m  odel  and  adjusting  the  quantity  of  elements  used  in  aver  aging  to 
determine  the  strain  at  a  gage  location  for  comparison  to  the  experimental  data. 
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VI.  CONCLUDING  REMARKS  AND  RECOMMENDATIONS 


Thin  com  posite  p  late  structures  were  exam  ined  under  various  con  ditions  to 
investigate  the  effect  of  FSI  on  dynam  ic  be  haviors.  Overall,  water  influenced 
significantly  both  kinetic  and  st  rain  energies  of  the  com  posite  structures  by  greatly 
reducing  their  m  agnitudes  and  frequencies.  The  FSI  greatly  suppressed  the  oscillatory 
nature  of  dynamic  responses  of  the  structures  .  Whether  a  structure  is  wet  on  one-side  or 
two-sides,  the  FSI  effect  was  very  clear  even  though  the  two-sides  wet  structures  showed 
a  greater  FSI  effect. 

The  boundary  condition,  either  clamped  or  simple,  has  sim  ilar  behaviors  and  is 
thereby  not  a  significant  contri  butor  for  FSI.  The  size  and  sh  ape  of  the  composite  plates 
was  shown  to  have  m  inor  differences  in  FSI.  The  m  ethod  of  loading  the  plate,  either 
concentrated  force,  uniform  pressure  or  impact,  showed  some  difference  on  the  degree  of 
FSI.  Interestingly,  the  shape  of  t  he  im  pacting  object  (contact  shape)  gave  different 
degrees  of  FSI  for  equivalent  im  pact  velocities.  The  larg  est  variation  of  FSI  was  due  to 
differences  in  m  aterial  properties  such  as  de  nsity  and  elastic  m  odulus.  As  a  result,  it  is 
critical  to  understand  and  incorporate  the  FSI  effects  when  designing  reliable  composite 
structures  employed  in  an  underwater  environment. 

Future  work  should  exam  ine  the  dynam  ic  behavior  of  composites  which  include 
moisture  ab  sorption  ef  fects.  Additiona  lly,  various  types  of  co  mposites  should  be 
compared  for  determ  ination  of  th  e  best  response  behavior  properties  and  m  inimal 
moisture  absorption.  Finally,  both  num  erical  and  experim  ental  work  should  be 
conducted  to  m  onitor  com  posite  behavior  in  failure.  The  failure  modes  should  be 
investigated  as  to  wheth  er  they  are  m  atrix  or  fiber  failure,  delam  ination  or  a  m  ixture  of 
failure  modes. 
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Norm.  Displacement,  Cntr  Top  Plate,  IWet  Norm.  Displacement,  Cntr  Top  Plate,  2Wet  Norm.  Displacement,  Cntr  Top  Plate,  Dry 


APPENDIX  A:  ADDITIONAL  FIGURES  FOR  CLAMPED  AND 
SIMPLE  BOUNDARY  WITH  CONCENTRATED  FORCE  LOAD 


Time  (sec.) 


Time  (sec.) 


Figure  55.  Displacement  and  Strain  Energy  Comparison  of  Clamped  and  Simple 

Boundary  with  Concentrated  Force  Load 


57 


Figure  56.  Normal  and  Shear  Strain  Comparison  at  Center  Position  for  Clamped 
versus  Simple  Boundary  with  Concentrated  Force  Load 
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Figure  57.  Normal  and  Shear  Strain  Comparison  at  Side  Position  for  Clamped  versus 

Simple  Boundary  with  Concentrated  Force  Load 
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APPENDIX  B:  ADDITIONAL  FIGURES  FOR  FORCE  AND 
PRESSURE  LOAD  COMPARISON  WITH  CLAMPED  BOUNDARY 


Figure  58.  Comparison  of  Dry  Structure  Response  for  Displacement,  Strain  and 
Kinetic  Energies  Between  Force  and  Pressure  Loading  with  Clamped  Boundary 


61 


Norm,  g^ ,  Qrtr  Bottom  Plate,  Dry  Norm,  g^ ,  Side  Bottom  Plate,  Dry  Norm,  g^ ,  Cntr  Bottom  Plate,  Dry 


Figure  59.  Normal  and  Shear  Strains  for  Comparison  of  Dry  Structure  with  Clamped 

Boundary  between  Force  and  Pressure  Loading 
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Norm,  e  ,  Qrtr  Bottom  Plate,  2Wet  Norm,  g  .Side  Bottom  Plate,  2Wet  Norm,  g  ,  Cntr  Bottom  Plate,  2Wet 
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Figure  60. 


Normal  and  Shear  Strains  for  Comparison  of  Two-sides  Wet  Structure 
with  Clamped  Boundary  between  Force  and  Pressure  Loading 
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Figure  61.  Normal  and  Shear  Strains  for  Comparison  of  One-side  Wet  Structure  with 
Clamped  Boundary  between  Force  and  Pressure  Loading 
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APPENDIX  C:  ADDITIONAL  FIGURES  FOR  FORCE  AND 
PRESSURE  LOAD  COMPARISON  WITH  SIMPLE  BOUNDARY 


Figure  62.  Comparison  of  Dry  Structure  Response  for  Displacement,  Strain  and 
Kinetic  Energies  between  Force  and  Pressure  Loading  with  Simple  Boundary 
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Figure  63.  Normal  and  Shear  Strains  for  Comparison  of  Dry  Structure  with  Simple 

Boundary  between  Force  and  Pressure  Loading 


66 


Norm,  e  ,  Qrtr  Bottom  Plate,  2Wet  Norm,  e  ,  Side  Bottom  Plate,  2Wet  Norm,  s  ,  Cntr  Bottom  Plate,  2Wet 


Figure  64. 


Normal  and  Shear  Strains  for  Comparison  of  Two-sides  Wet  Structure 
with  Simple  Boundary  between  Force  and  Pressure  Loading 
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Norm.g  ,  Qrtr  Bottom  Plate,  1  Wet  Norm.g  ,  Side  Bottom  Plate,  1  Wet  Norm.g  ,  Cntr  Bottom  Plate,  1  Wet 


Figure  65.  Normal  and  Shear  Strains  for  Comparison  of  One-side  Wet  Structure  with 
Simple  Boundary  between  Force  and  Pressure  Loading 
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APPENDIX  D:  ADDITIONAL  FIGURES  FOR  PLATE  SIZE 
EFFECTS  WITH  CONCENTRATED  FORCE  LOAD  AND 
CLAMPED  BOUNDARY 


Figure  66.  Comparison  of  Displacement  Response  for  Three  Structures  Due  to  Size 
Variation  Effects  with  Concentrated  Force  and  Clamped  Boundary 
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Figure  67.  Normal  and  Shear  Strains  for  Comparison  of  Differ  Plate  Sizes  for  Dry 
Structure  with  Concentrated  Force  and  Clamped  Boundary 
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Norm.  £o ,  Qrtr  Bottom  Plate,  2Wet  Norm,  ,  Side  Bottom  Plate,  2Wet  Norm,  e  Cntr  Bottom  Plate,  2Wet 


Figure  68.  Normal  and  Shear  Strains  for  Comparison  of  Differ  Plate  Sizes  for  Two- 
sides  Wet  Structure  with  Concentrated  Force  and  Clamped  Boundary 
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Norm.  s^.Qrtr  Bottom  Plate,  1  Wet  Norm. ,  Side  Bottom  Plate,  1  Wet  Norms  ,  Cntr  Bottom  Plate,  IWet 


Figure  69.  Normal  and  Shear  Strains  for  Comparison  of  Differ  Plate  Sizes  for  One- 
side  Wet  Structure  with  Concentrated  Force  and  Clamped  Boundary 
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APPENDIX  E:  ADDITIONAL  FIGURES  FOR  PLATE  SHAPE 
EFFECTS  WITH  CONCENTRATED  FORCE  LOAD  AND 
CLAMPED  BOUNDARY 


Figure  70.  Comparison  of  Displacement  Response  for  Three  Structures  Due  to  Shape 
Effects  with  Concentrated  Force  and  Clamped  Boundary 
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Figure  7 1 . 


Normal  and  Shear  Strains  for  Comparison  of  Differ  Plate  Shapes  for  Dry 
Structure  with  Concentrated  Force  and  Clamped  Boundary 


74 


Norm.  e^.Qrtr  Bottom  Plate,  2Wet  Norm.  E)(x ,  Side  Bottom  Plate,  2Wet  Norm.e  ,  Cntr  Bottom  Plate,  2Wet 


Figure  72.  Normal  and  Shear  Strains  for  Comparison  of  Differ  Plate  Shapes  for  Two- 
sides  Wet  Structure  with  Concentrated  Force  and  Clamped  Boundary 
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Norm. Qrtr Bottom  Plate,  1  Wet  Norm. g^,  Side  Bottom  Plate,  1  Wet  Norm.g  ,  Cntr Bottom  Plate,  1  Wet 


Figure  73. 


Normal  and  Shear  Strains  for  Comparison  of  Differ  Plate  Shapes  for  One- 
side  Wet  Structure  with  Concentrated  Force  and  Clamped  Boundary 


76 


APPENDIX  F:  ADDITIONAL  FIGURES  FOR  COMPOSITE 
DENSITY  EFFECTS  WITH  CONCENTRATED  FORCE  LOAD  AND 

CLAMPED  BOUNDARY 


The  following  use  composite  density  of  2020  kg/m  3  and  modulus  1.7e  10  GPa  for 
nonnalization  in  each  of  the  three  structures. 


a)  Dry  Structure 


Figure  74.  Comparison  of  Displacement  Response  for  Three  Structures  Due  to 
Density  Effects  with  Concentrated  Force  and  Clamped  Boundary 
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Figure  75.  Normal  and  Shear  Strains  for  Comparison  of  Differ  Density  for  Dry 
Structure  with  Concentrated  Force  and  Clamped  Boundary 
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Norm.  Exx  ,Qrtr  Bot.  Plate,  2Wet  Norm.  Exx  .Side  Bot.  Plate,  2Wet  Norm.E  ,Cntr  Bot.  Plate,  2Wet 


Time  (sec.) 


Time  (sec.) 


Figure  76.  Normal  and  Shear  Strains  for  Comparison  of  Differ  Density  for  Two-sides 
Wet  Structure  with  Concentrated  Force  and  Clamped  Boundary 
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Norm.  E)<x,QrtrBot.  Plate,  1  Wet  Norm.  E><x  .Side  Bot.  Plate,  1  Wet  Norm.E  ,Cntr  Bot.  Plate,  1  Wet 


Figure  77.  Normal  and  Shear  Strains  for  Comparison  of  Differ  Density  for  One-side 
Wet  Structure  with  Concentrated  Force  and  Clamped  Boundary 
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APPENDIX  G:  ADDITIONAL  FIGURES  FOR  COMPOSITE 
ELASTIC  MODULUS  EFFECTS  WITH  CONCENTRATED  FORCE 
LOAD  AND  CLAMPED  BOUNDARY 


The  following  use  composite  density  of  2020  kg/m  3  and  modulus  1.7e  10  GPa  for 
nonnalization  in  each  of  the  three  structures. 


a)  Dry  Structure 


b)  Two-sides  wet  Structure 


c)  One-side  Wet  Structure 


Figure  78.  Comparison  of  Displacement  Response  for  Three  Structures  Due  to 
Elastic  Modulus  Effects  with  Concentrated  Force  and  Clamped  Boundary 
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Figure  79.  Normal  and  Shear  Strains  for  Comparison  of  Differ  Elastic  Modulus  for 
Dry  Structure  with  Concentrated  Force  and  Clamped  Boundary 
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Norm.  Exx  ,Qrtr  Bot.  Plate,  2Wet  Norm.  E>s,  .Side  Bot.  Plate,  2Wet  Norm.  E  ]Cntr  Bot.  Plate,  2Wet 
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Figure  80.  Normal  and  Shear  Strains  for  Comparison  of  Differ  Elastic  Modulus  for 
Two-sides  Wet  Structure  with  Concentrated  Force  and  Clamped  Boundary 
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Norm.  e^.QrtrBot.  Plate,  1  Wet  Norm.  E><><  .Side  Bot.  Plate,  1  Wet  Norm.e  .CntrBot.  Plate,  1  Wet 


Figure  81.  Normal  and  Shear  Strains  for  Comparison  of  Differ  Elastic  Modulus  for 
One-side  Wet  Structure  with  Concentrated  Force  and  Clamped  Boundary 
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APPENDIX  H:  ADDITIONAL  FIGURES  FOR  IMP  ACTOR  SHAPRE 
EFFECTS  WITH  CLAMPED  BOUNDARY 


The  following  com  pare  circular  face  to  s  quare  face  im  pactor,  with  equal  im  pact 
area  and  equal  m  ass,  for  two-sides  and  one-sid  e  wet  structures  norm  alized  to  respectiv  e 
dry  structure. 


Two-sides  Wet  One-side  Wet 


Figure  82.  Comparison  of  Displacement  Response  for  Two-sides  and  One-side  Wet 

Structures  Due  to  Impactor  Shape  Effects 
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Figure  83.  Comparison  of  Strain  Energy  Response  for  Two-sides  and  One-side  Wet 

Structures  Due  to  Impactor  Shape  Effects 
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Figure  84.  Comparison  of  Kinetic  Energy  Response  for  Two-sides  and  One-side  Wet 

Structures  Due  to  Impactor  Shape  Effects 
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Figure  85.  Normal  and  Shear  Strain  Comparison  of  Different  Impactor  Shape  for 
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Two-sides  and  One-side  Wet  Structure  at  Center  Location 
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Figure  86.  Normal  and  Shear  Strain  Comparison  of  Different  Impactor  Shape  for 
Two-sides  and  One-side  Wet  Structure  at  Side  Location 
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Two-Sides  Wet 


One-side  Wet 


Figure  87.  Normal  and  Shear  Strain  Comparison  of  Different  Impactor  Shape  for 
Two-sides  and  One-side  Wet  Structure  at  Quarter  Location 
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APPENDIX  I:  ADDITIONAL  FIGURES  FOR  IMPACT  VELOCITY 

AND  SHAPE  EFFECTS 


The  following  com  pare  circular  face  to  s  quare  face  im  pactor,  with  equal  im  pact 
area  and  equal  mass,  for  different  velocities  of  the  three  structures,  normalized  to  1  m/s. 


Dry 


Time  (sec.) 


Circular  Square 

Comparison  of  Displacement  Response  for  Three  Structures  Due  to 
Impact  Velocity  Effects  for  Circular  and  Square  Faced  Impactor 
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Figure  88. 
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Figure  89.  Comparison  of  Strain  Energy  Response  for  Three  Structures  Due  to 
Impact  Velocity  Effects  for  Circular  and  Square  Faced  Impactor 
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Figure  90.  Comparison  of  Kinetic  Energy  Response  for  Three  Structures  Due  to 
Impact  Velocity  Effects  for  Circular  and  Square  Faced  Impactor 
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Figure  91.  Comparison  of  Normal  Strain  at  Center  Location  for  Three  Structures  Due 
to  Impact  Velocity  Effects  for  Circular  and  Square  Faced  Impactor 
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Figure  92.  Comparison  of  Shear  Strain  at  Center  Location  for  Three  Structures  Due 
to  Impact  Velocity  Effects  for  Circular  and  Square  Faced  Impactor 
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